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, , Abstract 

• Using a nonlocal extension of the SU(3) Nambu-Jona Lasinio model, which reproduces several of 

! the key features of Quantum Chromodynamics, we show that mixed phases of deconfined quarks 

I ' 

Q ' and confined hadrons (quark-hybrid matter) may exist in the cores of neutron stars as massive as 

in' 

c/3 ■ around 2.1 Mq. The radii of these objects are found to be in the canonical range of ~ 12 — 13 km. 

According to our study, the transition to pure quark matter does not occur in stable neutron stars, 
^ \ but is shifted to neutron stars which are unstable against radial oscillations. The implications of 

our study for the recently discovered, massive neutron star PSR J1614-2230, whose gravitational 
. mass is 1.97 it 0.04 M©, are that this neutron star may contain an extended region of quark-hybrid 



matter at it center, but no pure quark matter. 
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Introduction - Quantum Chromodynamics (QCD) has the properties of asymptotic free- 
dom and confinement. The former imphes that in the high- momentum transfer regime, the 
quarks behave essentially as free particles, i.e., the interaction between two quarks due to 
gluon exchange is very weak. This regime can therefore be treated using perturbation the- 
ory, where the quark-gluon coupling constant serves as an expansion parameter. For this 
momentum range, the dispersion processes can be calculated very accurately. By contrast, 
at low-momentum transfers (< 1 GeV) QCD becomes highly nonlinear, which prevents the 
use of perturbative methods. One of the renowned effective models that serves as a suit- 
able approximation to QCD in the low-energy regime is the quark version of the Nambu 
Jona-Lasinio (NJL) model l|-l3|. In this model chiral symmetry constraints are taken into 
account via effective interactions between quarks, through local four-point vertex interac- 
tions. The drawbacks of using local interactions are that the model must be regularized to 
avoid divergences in the loop integrals, and that the model is non- confining. The absence 
of confinement is essentially related to the fact that the dynamically generated constituent 
quark masses are momentum independent. Since the 1990 's, there have been investigations 
proposing nonlocal interactions to solve these problems ^. One interesting suggestion arises 
from the relationship between the NJL model and the model of one-gluon exchange where 
an effective gluon propagator is used to generate effective interactions between quarks. This 
provides a natural way to introduce a nonlocality in the quark-quark interaction, which can 
be characterized by a model-dependent form factor, g{p) j^. 

In this paper we analyze the global structure and composition of massive neutron stars 
in the framework of an extended version of the nonlocal SU(3) NJL model. Of particular 
interest is the question as to whether or not massive neutron stars, such as the recently 
discovered pulsar PSR J1614-2230 whose gravitational mass was found to be L97±0.04Mq 
6], may contain stellar cores made of deconfined quark matter tIB. [igI]. Tendentially, one 
could argue that quark deconfinement may not occur in the cores of such high-mass neutron 
stars, since the underlying nuclear equation of state must be extremely stiff in order to 
support such high mass neutron stars. As will be shown in this paper, arguments along 
these lines appear premature. 

Our study is based on the latest set of model parameters of the nonlocal SU(3) NJL 
model. The parameter fit is performed for a phenomenological value of the strange quark 
mass of TTis = 140.7 MeV. Details of the NJL formalism at zero chemical potential and finite 



temperature can be found in 10|-|l2[|. 

Over the last two decades, several authors have started to take into account the effect 
of .he ,ua.l. vecto. interaction in effective chiral .ode.s i.Ue NJlB g 
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-llTI]. It is known 



that the repulsive character of the vector coupling in these models affects the quark-hadron 
phase transition and moves the chiral restoration to a larger value of the quark chemical 
potential 15|. Thus, if the quark deconfined transition in the cores of neutron stars is 
modeled by a NJL-like model, it is expected that the vector coupling contribution modifies 
the nuclear equation of state (EoS) and hence the mass-radius relationship of neutron stars. 
Most of the NJL studies of neutron stars are treating the interactions among quarks in 
terms of local fermion-fermion couplings and/or impose the condition of local electric charge 



neutrality on the stellar matter js, 9|, ll6|. In this paper, we are using a generalized version 
of the NJL model where the interactions are nonlocal and momentum-dependent, and the 
condition of local charge neutrality is replaced with the more relaxed condition of global 
charge conservation, as neutron stars are characterized by two rather than one conserved 



charge 



18| . As a consequence, the pressure in the mixed quark-hadron phase varies with 



density and is, therefore, not a priori excluded from neutron stars. 



19|. However, we shall 



To include the vector interaction in the NJL model we follow 
consider three different vector fields, one for each quark flavor, instead of a single vector 
field for all quarks. 



Description of quark matter phase in the framework of the nonlocal SU(3) NJL model - We 
start from the Euclidean effective action associated with the nonlocal SU(3) quark model. 



d^x {ipix) [—i'j^dfj, + rh] tlj{x) 
-Y[^'{^^^'{x)+J^{x)j^{x)] 
-f T^^c K(x)jf (x)jf (x) - 3 fAx)j[{x)j^{x)] 



(1) 

where ip is a. chiral U{3) vector that includes the light quark fields, ip = [u d s)^, and 
m = diag(mu, m^, m^) stands for the current quark mass matrix. For simplicity we consider 
the isospin symmetry limit, in which m„ = = m. The fermion kinetic term includes the 
convariant derivative = d^ — lA^, where are color gauge fields, and the operator 7/^9^ 



in Euclidean space is defined as 7- V + 74^, with 74 = i'jQ. The currents jf'^(x) and jyj{x) 
are given by 



d'^z g{z) i){x + ^ \ai){x-^ , 
d^z g{z) ^ (a; + I) i IbK ^ " 



(x) = I d z g{z) i'f[x + -) Y i^f(x 



(2) 



2/ ' V 2 

where g{z) is a form factor responsible for the non-local character of the interaction, and 
the matrices A^, with a = 0, .., 8, are the usual eight Gell-Mann 3x3 matrices - generators 
of SU(3) - plus Aq = a/2/3 I3X3. Finally, the constants Tabc in the t'Hooft term accounting 
for flavor-mixing are defined by 



Tabc — 21 ^ijk ^mnl {^a)im i^b) jn i^c) f^i ■ 



(3) 



After standard bosonization of Eq. ([T]), the integrals over the quark fields can be performed 
in the framework of the Euclidean four-momentum formalism. The grand canonical potential 
in the mean-field approximation at zero temperature, including the vector coupling, is then 
given by 
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where Nc = 3, Ef = \ p"^ + rn^r and we have defined 



f = {po + + r 

The masses of free quarks are denoted by m/, where / 
related to the vector current density jyj of Eq. ([2]). 



(5) 

u, d, s. The wyj mean field is 



The momentum-dependent constituent quark masses Mj depend explicitly on the quark 
mean fields d"/, 

Mfiuj) = mf + afg{uj), (6) 
where gioo"^) denotes the Fourier transform of the form factor g{z). 



Following Ref. 19|, the quark vector interaction shifts the quark chemical potential ac- 
cording to 

Jlf = fif- wvj. (7) 

Note that the shifting of the quark chemical potential does not affect the nonlocal form 
factor g{ujj) as discussed in 19|. The mean-field values of the auxiliary fields Sf are given 
by 



io| 
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(8) 



In this paper we adopt a Gaussian form for the nonlocal form factor g, 



g{u') = exp i-oo'/A') 



(9) 



where A plays a role for the stiffness of the chiral transition. This parameter, together with 
the current quark mass m of up and down quarks and the coupling constants Gg and H in 
Eq. (jlj), have been fitted to the pion decay constant, and meson masses m^, m^, and m^/, 
1^. The resuh of this fit is m = 6.2 MeV, A = 706.0 MeV, G.A^ = 15.04, 
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as described in 

ifA^ = —337.71. The strange quark current mass is treated as a free parameter and was 
set to rris = 140.7 MeV. The strength of the vector interaction Gy is usually expressed in 
terms of the strong coupling constant Gs- To account for the uncertainty in the theoretical 
predictions for the ratio Gy/Gs, we treat the vector coupling constant as a free parameter 



20N22| . which varies from to 0.1 Gg- 

Using these parametrizations, the fields cx/ and ^v,f can be determined by minimizing 



Eq. 
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(10) 



Description of confined hadronic matter - The hadronic phase is described in the framework 
of the non-linear relativistic mean field theory 



23 



26| . where baryons (neutrons, protons. 



hyperons) interact via the exchange of scalar, vector and isovector mesons (cr, cj, p, respec- 



tively). The Lagrangian of the theory is given by 

£ = 



S=n,p,A,E,H 



- {rriN - gaO-)\ilJB + - {d^adf'a - m^a ) 



1 



(11) 



III/ 



\=e~ 



where B sums all baryon states which become populated in neutron star matter 25|, l26 |. 
The quantities gp, g^j, and g^^ are the meson-baryon coupling constants. Non-linear a-meson 
-interactions are taken into account in Eq. fill I) via the terms proportional to and c^r 
261 ] ■ We have solved the equations of motion for the baryon and meson field equations, 



sel 



25 



which follow from Eq. ffTTl) . for the relativistic mean-field approximation 25|, |26|. For this 
approximation the meson fields cr, u, p are approximated by their respective mean- field 
values o = (cr), Uj = (u), and p = (poa) 25|, l26|]. The parameters of the model, labeled GMl, 



are adjusted to the properties of nuclear matter at saturation density. They are taken from 

The condition of weak equilibrium requires the presence of electrons and muons, which 
are treated as free relativistic quantum gases, as described by the last term on the right- 
hand-side of Eq. ( ITT]) . Neutron star matter is characterized by the conservation of electric 
and baryon number. This feature leads to the chemical equilibrium condition 



pi -Bi Pn Qi pe ) 



(12) 



where pn and pe denote the chemical potentials of neutrons and electrons, respectively. The 
quantities Bi and Qi stand for the baryon numbers and the electric charges of the mesons 
and baryons of Eq. (11 ip . Equation (11 2p greatly simplifies the mathematical analysis, since 
only the knowledge of two independent chemical potentials, pn and pe, is necessary. The 



latter are obtained from 



25 



26| 
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FIG. 1. (Color online) Pressure, P (solid lines), baryon chemical potential, (dashed lines), 
and electron chemical potential, /Xg (dotted lines) as a function of baryon number density, p, in 
units of the normal nuclear matter density, pQ = 0.16 fm~^. The hatched areas denote the mixed 
phase regions where confined hadronic matter and deconfined quark matter coexist. The solid dots 
indicate the central densities of the associated maximum-mass stars, shown in Fig. [SI and x is the 
respective fraction of quark matter inside of them. The results are computed for three different 
values of the vector coupling constant, ranging from 0, to 0.05 Gg, to 0.1 Gs- 

where m*^ = — g^^a denote the effective medium- modified baryon masses, ks and kx are 
the Fermi momenta of baryons and leptons, respectively, and is the third component of the 
isospin vector of a baryon of type B. Finally, aside from chemical equilibrium, the condition 
of electric charge neutrality for confined hadronic matter is also of critical importance for 
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FIG. 2. (Color online) Volume fraction, X; of quark phase as a function of baryon number density, 
p, in units of normal nuclear matter density, po = 0.16 fm~^. The solid dots indicate the central 
densities of the respective maximum-mass stars shown in Fig. [31 



the composition of neutron star matter. This condition is given by 25|, |26| 



B A 

where Jb denotes the spin of baryon B. 

For the quark phase, the chemical potentials associated with quarks and electrons follow 
from Eq. ^ as fJ'u = f^b ~ 2/ie/3 and fi^ = f^s = f^b ~\~ /^e/3, where /if, = /i„/3 stands for the 
baryon chemical potential, expressed in terms of the chemical potential of neutrons. 

Description of the mixed phase of quarks and hadrons - To determine the mixed phase region 
of quarks and hadrons, we start from the Gibbs condition for phase equilibrium between 
hadronic (H) and quark (Q) matter, 

Pff /ie, {0}) = PgifJ-a, /^e) , (15) 

where P„ and P, denote the p.e=swes of had.cnc .natter and qua* matte, respectively 
[18( 1 . The quantity {0} in Eq. (HM stands collectively for the field variables (a, u, p) and 
Fermi momenta (/c^, kx) that characterize a solution to the equations of confined hadronic 
matter. We use the symbol x = ^q/^ to denote the volume proportion of quark matter, 
Vq, in the unknown volume V. By definition, x then varies between and 1, depending 
on how much confined hadronic matter has been converted to quark matter. Equation (fT5|) 
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is to be supplemented with the conditions of global baryon charge conservation and global 



electric charge conservation. The global conservation of baryon charge is expressed as 



18| 



Pb = X Pgif^n, Pe) + (1 - X) PH{Pn, Pe, {4>}) , (16) 

where pg and pn denote the baryon number densities of the quark phase and hadronic phase, 
respectively. The global neutrality of electric charge is given by 



= X (iQiPn, Pe) + (1 - X) quiPn, Pe, {0}) , (17) 

with qq and qn denoting the electric charge densities of the quark phase and hadronic phase, 
respectively. We have chosen global rather than local electric charge neutrality, since the 
latter is not fully consistent with the Einstein-Maxwell equations and the micro physical 
condition of /3-equilibrium and relativistic quantum statistics, as shown in 
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. Local NJL 



studies carried out for local electric charge neutrality have been reported recently in Refs. 
js, 9|. In Ref. jo] the nonlinear relativistic mean-field model GMl and the local SU(3) NJL 
model with vector interaction were used to describe the hadronic and quark matter phases, 
respectively. The authors found that the observation of neutron stars with masses a few 
percent higher than the 1.97 ± 0.04 M0 would be hard to explain unless, instead of GMl, 
one uses a very stiff model for the hadronic EOS with nucleons only. They obtain neutron 
stars with stable pure quark matter cores in their centers. In contrast to this, we find that 
such neutron stars are not be stable if the nonlocal NJL model is used instead of the local 
model and the less stringent condition of global electric charge neutrality is imposed on the 
composition of the stellar matter. 

The results for the mixed phase region are shown in Figs. [1] and |2j Our calculations show 
that the inclusion of the quark vector coupling contribution shifts the onset of the phase 
transition to higher densities, and also narrows the width of the mixed quark-hadron phase, 
when compared to the case Gy = 0. The mixed phases range from 3.2 — 8.2po, 3.8 — 8.5^0, 
and 4.5 — 8.9po for vector coupling constants Gy/Gg = 0, 0.05, 0.1, respectively. We note 



that there is considerable theoretical uncertainty in the ratio of Gy / Gg 29| since a rigorous 
derivation of the effective couplings from QCD is not possible. Combining the ratios of 
Gv/Gg from the molecular instanton liquid model and from the Fierz transformation, the 



value of Gy/Gs is expected to be in the range ^ Gy/Gg ^ 0.5 [30|. For our model, values 
of Gy/Gg > 0.1 shift the onset of the quark-hadron phase transition to such high densities 
that quark deconfinement can not longer occur in the cores of neutron stars. 
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Next we determine the bulk properties of spherically symmetric neutron stars for the 
collection of equations of state shown in Fig. [1] The properties are determined by solving 
the Tolmann-Oppenheimer-Volkoff (TOV) equation of general relativity theory 31|. The 
outcome is shown in Fig. [31 One sees that depending on the value of the vector coupling 
constant, Gy, the maximum neutron star masses increase from 1.87 Mq for Gy = 0, to 
2.00 Mq for Gv = 0.05 G„ to 2.07 Mq for Gy = 0.1 G,. The heavier stars of all three stellar 
sequences contain mixed phases of quarks and hadrons in their centers. The densities in 
these stars are however not high enough to generate pure quark matter in the cores. Such 
matter forms only in neutron stars which are already located on the gravitationally unstable 
branch of the mass-radius relationships. Another intriguing finding is that neutron stars 
with canonical masses of around 1.4 Mq do not posses a mixed phase of quarks and hadrons 
but are made entirely of confined hadronic matter. 



Summary and Conclusions - In this paper, we show that high-mass neutron star, such as 
PSR J1614-2230 with a gravitational mass of 1.97 ± 0.04 j^, may contain mixtures of 
quarks and hadrons in their central regions. Our analysis is based on a nonlocal extension 
of the SU(3) Nambu-Jona Lasinio model, which reproduces some of the key features of 
Quantum Chromodynamics at densities relevant to neutron stars. Critical is the inclusion 
of the quark vector coupling contribution in the nonlocal SU(3) NJL model. Our results also 
show that the transition to pure quark matter occurs only in neutron stars which lie already 
on the gravitationally unstable branch of the mass-radius relationship. The existence of 
pure quark matter in massive (as well as in all other, less massive) neutron stars would thus 
be ruled out by our study. 
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(a) Mass-radius relationships of neutron stars (b) Enlargement of the circled region of Fig. 3(a) 
with quark-hybrid matter in their centers. The labels 'MP' and 'QP' stand for mixed phase 

and pure quark phase, respectively. The solid 
vertical bars denote maximum-mass stars. 



FIG. 3. (Color online) Depending on the strength of the vector repulsion {Gy) of the nonlocal 
NJL model, maximum masses up to 2.1 Mq are obtained. With increasing stellar mass, the stellar 
core compositions consist of either only nucleons, nucleons and hyperons, a mixed phase of quarks 
and hadrons (MP), or a pure quark matter phase (QP). The latter, however, exists only in neutron 
stars which lie to the left of their respective mass peaks. Such stars are unstable against radial 
oscillations and thus cannot exist stably in the universe. In contrast to this, all neutron stars on 
the MP branches up to the mass peaks are stable. The gray band denotes the 1-a error bar of the 
M = (1.97 ± O.O4)M0 neutron star PSR J1614-2230 [e]. 
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